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The ubiquitous presence of dark matter in the universe is today a central tenet in modern 
cosmology and astrophysiciP. Ranging from the smallest galaxies to the observable universe, 
the evidence for dark matter is compelling in dwarfs, spiral galaxies, galaxy clusters as well as 
at cosmological scales. However, it has been historically difficult to pin down the dark matter 
contribution to the total mass density in the Milky Way, particularly in the innermost regions 
of the Galaxy and in the solar neighbourhoo(P. Here we present an up-to-date compilation 
of Milky Way rotation curve measurementiP®!, and compare it with state-of-the-art bary- 
onic mass distribution modeli^i^^. We show that current data strongly disfavour baryons as 
the sole contribution to the galactic mass budget, even inside the solar circle. Our findings 
demonstrate the existence of dark matter in the inner Galaxy while making no assumptions 
on its distribution. We anticipate that this result will compel new model-independent con¬ 
straints on the dark matter local density and profile, thus reducing uncertainties on direct 
and indirect dark matter searches, and will shed new light on the structure and evolution of 
the Galaxy. 

Existing studies of the dark matter density in the inner Galaxy fall into two eategories: mass¬ 
modelling and loeal measurements. In mass-modelling, the distribution of dark matter is assumed 
to follow a density profile inspired by numerieal simulations, typieally an analytie fit sueh as the 
well-known Navarro-Frenk-WhiteP^or EinastcP^ profiles, with two or more free parameters whose 
best-fit values are then determined from dynamieal eonstraints. The statistieal error on the dark 
matter density in the inner Galaxy - and in partieular in the solar neighbourhood - is in this ease 
very small, of order 10%P^, but this only refleets the strong assumptions made about the distri¬ 
bution of dark matter. The latter is in faet observationally unknown, and the aforementioned 
elasses of profiles are inspired by simulations without baryons, whose role is not negligible in 
the inner Galaxy. Local measurements are instead based on the study of observables in the solar 
neighbourhoocP. These methods ean be used to assess the evidenee for dark matter loeally through 
an estimate of the gravitational potential from the kinematies of stars. However, the value found 
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for the local dark matter density is usually compatible with zero at ~3 ct unless one makes strong 
assumptions about the dynamics of the tracer populations. 

Here we report on a comparison of the observed rotation curve of the Galaxy with that 
expected from visible matter only. As we shall see, this approach provides an alternative way 
to constrain additional contributions of matter to the rotation curve, and therefore to infer the 
existence and abundance of dark matter. Although this has been historically one of the first methods 
to detect dark matter in external galaxies, it has long been thought to provide weak constraints in 
the innermost regions of the Milky Way, due to a combination of poor rotation curve data and 
large uncertainties associated with the distribution of baryons. We show that recent improvements 
on both fronts permit to obtain a convincing proof of the existence of dark matter inside the solar 
circle. 

We start by presenting a new, comprehensive compilation of rotation curve data derived 
from kinematic tracers of the galactic potential, which considerably improves upon earlier (partial) 
compilationd^^EH. Optimised to galactocentric distances i? = 3 — 20 kpc, our database includes 
gas kinematics (HI terminal velocitied^ll, CO terminal velocities, HI thicknesS, HII regionS®, 
giant molecular cloudS), star ki nematics (open clusterP, planetary nebulaeP^, classical cepheidS^, 
carbon star^^ and maser^. This represents an exhaustive survey of the literature that intention¬ 
ally excludes objects with kinematic distances only, and those for which asymmetric drift or large 
random motions are relevant. In total we have compiled 2780 measurements, of which 2174, 506 
and 100 from gas kinematics, star ki nematics and masers, respectively (see supplementary infor¬ 
mation). For each measurement, we translate the kinematic data into a constraint on the angular 
velocity ujc = vjR and on the galactocentric distance R. The upper panel of Fig. 1 shows the 
rotation curve Vc{R) for the full compilation of data. 

The contribution of stars and gas to the total mass of the Galaxy has historically been subject 
to significant uncertainties, in particular towards the innermost regions where its dynamical contri¬ 
bution is most important. Substantial progress has been made recently, and data-based models that 
encode the three-dimensional morphology of the baryonic distribution have become available in 
the literature. In order to bracket the uncertainties on the stellar and gas distribution, we consider 
here all possible combinations of a set of detailed models for the stellar bulgef^^^, stellar diskf^^^ 
and ga^ ^°F* l(see supplementary information). The stellar bulge models encompass alternative den¬ 
sity profiles in the inner few kiloparsecs and different configurations of the galactic bar. The stellar 
disks implemented provide instead the best descriptions of star observations across the Galaxy, in¬ 
cluding parametrisations with and without separation into thin and thick populations. Finally, the 
gas is split into its molecular, atomic (cold, warm) and ionised (warm, hot, very hot) components, 
paying special attention to the localised features in the range i? = 10 pc—20 kpc. 

The gravitational potential of each model is computed through multipole expansiorP^ (see 
supplementary information), and the corresponding rotation curve is shown in the lower panel of 
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Fig. 1 with its original normalisation. We calibrated each bulge with the microlensing optical depth 
measurement towards the central galactic region (r) = 2.17lo;3g ^ each disk with the 

constraint on the local surface density S* = 38 ± 4 Mq/pc®, and for the g as we adopted a CO- 
to-H 2 conversion factor of (0.5 — 3.0) x 10^° cm“^(K km/s)“^ for i? > 2 kpcJ^^Oll. This procedure 
ensures that all baryonic models comply with the existing observational constraints and moreover 
it assigns a realistic uncertainty to the contribution of each model to the rotation curve. 

We assess the evidence for an unseen (dark) component of the gravitational potential of our 
Galaxy in the form of a discrepancy between the observed rotation curve and that expected from the 
set of baryonic models described above. We stress that we do not make any assumption about the 
nature or distribution of dark matter: our analysis therefore provides a model-independent estimate 
of the amount of dark matter in the Galaxy. For each baryonic model, the two-dimensional chi- 
squareP is computed and used to assess the goodness-of-fit. We have explicitly checked through 
Monte Carlo calculations that this statistic has an approximate distribution for the case at hand. 
The analysis is restricted to galactocentric distances R > Rent = 2.5 kpc, below which the orbits 
of the kinematic tracers are significantly non-circular. We adopt a distance to the galactic centre 
Rq = 8 kpc, a local circular velocity no = 230 km/s, and a peculiar solar motiorl^ {U, V, W)^ = 
(11.10,12.24,7.25) km/s. 

The upper panel of Fig. 2 shows the angular velocity as a function of the galactocentric dis¬ 
tance. Observational data are shown with red dots, while the grey band shows the envelope of all 
baryonic models discussed above, which we interpret here as bracketing the possible contribution 
of baryons to the rotation curve. The discrepancy between observations and the expected contribu¬ 
tion from baryons is evident along the whole range of galactocentric distances above 6 — 7 kpc. The 
residuals are plotted in the middle panel of Fig. 2 for a fiducial baryonic model (the one shown by 
the solid black line in the upper panel), and they can be readily interpreted as the contribution of an 
extra component to the Newtonian gravitational potential of our Galaxy. Interestingly, the gravita¬ 
tional potential from a dark matter distribution such as those suggested by numerical simulations 
(Navarro-Frenk-White or Einasto profiles) smoothly fills the gap without fine tuning. 

The main conclusion of our analysis is summarised in the bottom panel of Fig. 2, where we 
plot the per degree of freedom for each baryonic model and for all data up to a given distance 
R (but above Rent)- The evidence for a dark component rises above 5a (thick red line) well 
inside the solar circle for all baryonic models. Indeed, whereas the relative discrepancy between 
observational data and baryonic models is higher at larger galactocentric distances, it is at lower 
distances that uncertainties are smallest. Hence, the evidence grows swiftly at relatively small 
distances and saturates at larger distances. We have tested the robustness of our results against 
variations of Ro,vo, peculiar solar motion, binning and data selection as well as against systematics 
due to spiral armd^. The results change only mildly for all cases, and the conclusions drawn from 
Fig. 2 remain unchanged (see supplementary information). 
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The comparison of the Milky Way observed rotation curve with the predictions of a wide 
array of baryonic models points strongly to the existence of a contribution to the gravitational po¬ 
tential of the Galaxy from an unseen, diffuse component. The statistical evidence is very strong 
already at small galactocentric distances, and it is robust against uncertainties on galactic morphol¬ 
ogy and kinematics. Without any assumption about the nature of this dark component of matter, 
our results open a new avenue for the determination of its distribution inside the Galaxy. This has 
powerful implications both on studies aimed at understanding the structure and evolution of the 
Milky Way in a cosmological context, and on direct and indirect dark matter searches, aimed at 
understanding the very nature of dark matter. 
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Figure 1: The rotation curve of the Milky Way. In the top panel we show our compilation of 
rotation curve observations as a function of galactocentric distance, including data from gas kine¬ 
matics (blue dots; HI terminal velocities, CO terminal velocities, HI thickness, HII regions, giant 
molecular clouds), star ki nematics (open green squares; open clusters, planetary nebulae, classical 
cepheids, carbon stars) and masers (open black circles). Error bars correspond to la uncertainties. 
The bottom panel displays the contribution to the rotation curve as predicted from different mod¬ 
els for the stellar bulge (blue), stellar disk (green) and gas (black). We assume a distance to the 
galactic centre Rq = 8 kpc in both panels, and a local circular velocity vq = 230 km/s in the top 
panel. 
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Figure 2: Evidence for dark matter. In the top panel we show the angular velocity measurements 
from the compilation shown in Fig. 1 (red dots) together with the bracketing of the contribution 
of all baryonic models (grey band) as a function of galactocentric distance. Error bars correspond 
to Icr uncertainties, while the grey band shows the envelope of all baryonic models including la 
uncertainties. The contribution of a fiducial baryonic model is marked with the black line. The 
residuals between observed and predicted angular velocities for this baryonic model are shown in 
the middle panel. The dashed blue line shows the contribution of a Navarro-Frenk-White profile 
with scale radius of 20 kpc normalised to a local dark matter density of 0.4 GeV/cm^. The bottom 
panel displays the cumulative reduced for o^ch baryonic model as a function of galactocentric 
distance. The black line shows the case of the fiducial model plotted in black in the top panel, 
while the thick red line represents the reduced corresponding to 5a significance. In this figure 
we assume a distance to the galactic centre Rq = 8 kpc and a local circular velocity vq = 230 
km/s, and we ignore all measurements below Rent = 2.5 kpc. 
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Supplementary Information 


Materials and Methods 

Rotation curve data 


The main features of the eompilation of rotation eurve data used in this work are summarised 
in Table SI. Every adopted objeet or region has assoeiated galaetie eoordinates (£, b), helioeentrie 
distanee d and helioeentrie line-of-sight veloeity The latter is usually reported in a given loeal 
standard of rest (LSR) frame rather than in the helioeentrie frame, so that one has to subtraet the 
peeuliar solar motion used in the original reference to find and then apply the desired peculiar 
solar motion (t/, V, W)^ to get the LSR line-of-sight velocity vl°^. We follow elosely each source 
referenee to assign errors to d and vl°^ and to aeeount for any peculiar motion associated with 
speeific objeets. Uneertainties on i and b are largely sub-dominant in all cases and are therefore 
negleeted. Finally, we exelude objects with insufficient or defieient data (e.g. on the distanee 
determination), too elose to the direetion of the galaetie eentre or anti galaetie eentre and any other 
objeets elassified as suspeet in the original referenees. After this seleetion, the eompilation eonsists 
of 2780 individual measurements with the breakdown shown in Table SI. 

We then eonstrain the rotation eurve of our Galaxy Vc{R) for any given ehoiee of the distanee 
to the galaetie eentre Rq and loeal cireular veloeity vq = Vc{Ro). Assuming eireular orbits for the 
objeets observed (a reasonable approximation outside the influenee of the galaetie bulge, i.e. R > 
Rent = 2.5 kpc), 

v\Tr = cos6 sinf , (1) 

where R = {d? cos^ b + R^ — 2RQd cos b cos For the partieular ease of terminal veloeities, 
6 = 0 and R = Rq \ sin£ |. When proper motions are available (e.g. for open elusters and masers), 
similar expressions apply for the objeet’s veloeity along the longitude and latitude directions. For 
each object in the eompilation, we have a measurement of R and we invert Eq. Q to obtain the 
angular cireular velocity Wc{R) = Vc{R)/R and propagated error. Notiee that we make use of 
the angular cireular velocity Wc rather than the actual circular velocity Vc, because the error of the 
latter is strongly correlated with the error of R. Instead, the errors of Wc and R are uncorrelated. As 
notieed long ageP, using Vc would introduee unneeessary eomplieations in the statistieal analysis 
(see below) and lead to a degradation of the aeeuraey. 

Baryonie modelling 


The exaet distribution of baryons in our Galaxy is not preeisely determined as of today, mak¬ 
ing it a major souree of uneertainty in the present study. There are three main baryonie eompo- 
nents: stellar bulge, stellar disk and gas. We have surveyed the literature exhaustively and eolleeted 
a wide range of data-based, three-dimensional morphologies for eaeh eomponent. This allows for 
a quantitative assessment of the bracketing due to baryonie modelling, as shown in Fig. 2 in the 
main text. The details of bulge, disk and gas models are given below. 
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The inner few kpe of the Milky Way are dominated by a triaxial, bar-shaped bulge of 
stard ^"^l^^ l ^ . Observations elearly plaee the near end of the bar at positive galaetie longitudes, 
but its preeise orientation and morphology are less eertain. For instanee, the distribution of red 
elump giants in the bulgeP^ is well fitted either by exponential or gaussian profiles (so-ealled E2 
and G2 models, respeetively). Apart from these two eonfigurations, we also eonsider alternative 
truneated power-law bulge^^lllll and a bar-shaped model with a nuelear eomponenP^. In view of 
reeent developments, the possibility of an extra (long) bai^and a double-ellipsoid bulgeP^are im¬ 
plemented as well. The normalisations of all seven models (and eorresponding uneertainties) are 
fixed by matehing the predieted mierolensing optieal depth towards (f, b) = (1.50°, —2.68°) to the 
2005 MACHO measurement (r) = 2.17^g gg x 10“^23l2l. The mierolensing eontribution due to 
disk stars is self-eonsistently aeeounted for in aeeordanee to the disk models deseribed below. 

The stellar disk has been modelled by different authors with the help of eomprehensive sur¬ 
veys of photometrie data aeross the Galaxy. Typieal parameterisations inelude thin and thiek disk 
populations, usually featuring double-exponential profiles. We eonsider alternative thin plus thiek 
eonfigurationd^31II as well as eonfigurations with a stellar halo eomponenfi^ni. The single max¬ 
imal disk suggested reeentl)^^^ is also implemented. All five models are normalised to the latest 
loeal surfaee density eonstraint S* = 38 ± 4 MQ/pe^, from whieh we propagate the uneertainty 
to the disk eomponent. 

Finally, a non-negligible part of the baryons in the Milky Way is in the form of gas, namely 
moleeular, atomie and ionised hydrogen and heavier elements. The distribution of eaeh eomponent 
is relatively well-known but extremely irregular. This is for example the ease of the gas within 
10 pe of the galaetie eentreP^, whieh for our purposes ean be safely eonsidered as a point-like 
mass. For the inner 2 kpe, we model moleeular and atomie hydrogen in the eentral moleeular zone 
and holed disk, and the distribution of ionised hydrogen is split into its warm, hot and very hot 
phased^. In the range i? = 2 — 20 kpe, instead, two alternative morphologied22EI] are used for 
eaeh gas eomponent. We set up in this way our two gas models, whose uneertainties are assigned 
by taking a CO-to-H 2 eonversion faetor of (2.5 — 10) x 10^® om“^(K km/s)“^ for R < 2 kpe and 
(0.5 — 3.0) X 10^° em“^(Kkm/s)“^ for i? > 2 kpd^^lkl. 

Onee a model for bulge, disk and gas is speeified, the individual gravitational potentials (and 
thus the individual eontributions to the rotation eurve) ean be easily eomputed through multipole 
expansiorP^. Expanding up to Imax = 2 (see below for a eonvergenee test) and averaging over 
the azimuthal direetion, we ean then derive the overall baryonie eontribution to the rotation eurve 

^ bulge + ^disk + ^las ^^e eorrespoudiug propagated uneertainties. 

Statistieal analysis 


The eentral task in this work is to eompare the observed rotation eurve of the Galaxy uJc{R) to 
that expeeted from baryons uJb{R), and deeide whether is eompatible with zero or not. This 

task is eomplieated by the sizeable error on R in the rotation eurve data, espeeially at intermediate 
and large R as shown in Eig. 1 of the main text. Given the large amount of observations (see 
Table SI), a eustomary teehnique usually adopted in the literature is that of binning the data. Sinee 
binning entails loss of information, we opt not to do it and use instead the full power of the data 
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taking proper account of galactocentric distance errors. (We did however check explicitly that a 
binned analysis with a weighted mean of the measurements just reinforces the results presented in 
the main text.) Following Ref. and introducing the reduced variables x = R/RodMAy = u/uq — I 
(with uq = vq/Ro), the two-dimensional reads 
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where {xb,i, yb,i) is the point in the baryonic curve yb{x) = uJb{R = xRq)/ujq — 1 that minimises 
di- Notice that the expression above can be applied because the errors on R and oj (i.e. x and 
y) are uncorrelated and hence the error ellipse is not tilted in the (i?, oj) plane. This would not 
be the case for the error-correlated pair {R,Vc). We have performed Monte Carlo calculations 
for a fiducial baryonic model and typical uncertainties ax,i, CTy^i, and verified that the statistic in 
Eq. Q follows approximately a distribution. The lower panel in Fig. 2 in the main text shows 
the reduced chi-square where the sum in Eq. (|^ is restricted to objects with galactocentric 

distances below a given R. Due to the breakdown of the assumption of circular orbits, all objects 
with R < Rent = 2.5 kpc are ignored in the analysis. 


Robustness of the results 


Our main findings are presented in Fig. 2 for a wide bracketing of baryonic models. This 
illustrates already the robustness of the results. However, Fig. 2 was obtained for fixed (albeit 
reasonable) choices of galactic parameters and data selection. Here we show that different choices 
lead to the same conclusion as in the main text, i.e. that baryons cannot explain alone the observed 
rotation curve in the inner Galaxy. 

We start by varying the fundamental galactic parameters, namely Rq, vq and {U^V^W)^. 
The existing determinations of Rq (e.g.l ^^F^l Ell) do not all agree with each other, but are confined to 
a reasonably narrow range, 8.0 ± 0.25 kpeP^. We a dopt a more conservative range 8.0 ± 0.5 kpc in 
our tests. A similar situation holds for vq (e.g.E3511EZl) where the range 230 ± 20 km/s encompasses 
most determinations. As for the peculiar solar motion, we take ((7, V, W)^ = (11.10,12.2 4, 7.2 5) 
km/spfl as our fiducial values, but consider as well other determinations in the literatureC^MZMI: 
{U,V,W)q = (10.00,5.25,7.17) km/s, ([/, E, 1E)q = (10.00,26.00,7.25) km/s (with uq = 218 
km/s) and (f/, V, W)^ = (10.7,15.6, 8.9) km/s (with vq = 240 km/s). The effect of taking different 
choices of galactic parameters is shown in Fig. SI. Also shown is the evidence obtained when 
using gas kinematics, star kinematics and masers separately. Furthermore, the convergence of the 
multipole expansion used to compute ujb was tested up to imax = 8. Finally, we checked the impact 
of streaming motions due to spiral arms by adding a 11.8 km/s systematic^ to each rotation curve 
measurement. In all cases the discrepancy between observed and predicted rotation curves is high 
(> 5 ct) and very robust against galactic nuisances, data selection and systematics. 
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Figure S1: The evidenee for dark matter in the inner Galaxy and its dependenee on galaetie param¬ 
eters, data seleetion and systematies. The plots show the eumulative redueed for eaeh baryonie 
model as a funetion of galaetoeentrie distanee. The thiek red line represents the redueed corre¬ 
sponding to 5cr significance, while the black line shows the case of the fiducial baryonie model. 
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Object type 

R [kpc] 

quadrants 

# objects 

HI terminal velocities 

Fich+ ’8P 

1 

OO 

b 

1,4 

149 

Malhotra ’9^ 

2.1 -7.5 

1,4 

110 

McClure-Griffiths & Dickey ’0"^ 

2.8-7.6 

4 

701 

HI thickness method 

Honma & Sofue ’9"^ 

6.8-20.2 

— 

13 

CO terminal velocities 

Burton & Gordon ’7d^ 

1.4-7.9 

1 

284 

Clemens ’8^ 

b 

1 

OO 

b 

1 

143 

Knapp+ ’8^ 

p 

a\ 

1 

bo 

1 

37 

Luna+ 

2.0-8.0 

4 

272 

HII regions 

Blitz ’7923 

8.7-11.0 

2,3 

3 

Fich+ ’8P 

9.4-12.5 

3 

5 

Turbide & Moffat ’93l23 

11.8-14.7 

3 

5 

Brand & Blitz ’93pl 

5.2-16.5 

1,2,3,4 

148 

Hou+ ’0#l 

3.5-15.5 

1,2,3,4 

274 

giant molecular clouds 

Hou+ ’O# 

6.0-13.7 

1,2,3,4 

30 

open clusters 

Frinchaboy & Majewski ’Oo^ 

4.6-10.7 

1,2,3,4 

60 

planetary nebulae 

Durand+ ’98121 

3.6-12.6 

1,2,3,4 

79 

classical cepheids 

Pont+ ’9423 

5.1 - 14.4 

1,2,3,4 

245 

Pont+ ’9'P 

10.2-18.5 

2,3,4 

32 

carbon stars 

Demers & Battinelli ’0"^ 

9.3-22.2 

1,2,3 

55 

Battinelli+ ’13P3 

12.1-24.8 

1,2 

35 

masers 

Reid+ ’ 1^ 

4.0-15.6 

1,2,3,4 

80 

Honma+ ’ li23 

7.7-9.9 

1,2,3,4 

11 

Stepanishchev & Bobylev ’ 1 1^ 

8.3 

3 

1 

Xu+ ’IsP 

7.9 

4 

1 

Bobylev & Bajkova ’Isl^Sl 

4.7 - 9.4 

1,2,4 
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Table SI: Our compilation of rotation curve data for the Milky Way. For each object type and 
reference, we report the range of galactocentric distance R (assuming a distance to the galactic 
centre Rq = 8 kpc), the galactic quadrant(s) as well as the number of objects analysed (after cuts). 
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